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The microscopic description of neutron scattering by 16 below 30 MeV is carried out by means of the 
continuum particle-vibration coupling (cPVC) method with the Skyrme nucleon-nucleon (NN) effective inter- 
action. In the cPVC method, a proper boundary condition on a nucleon in continuum states is imposed, which 
enables one to evaluate the transition matrix in a straightforward manner. Experimental data of the total and 
total-elastic cross sections are reproduced quite well by the cPVC method. An important feature of the result 
is the fragmentation of the single-particle resonance into many peaks as well as the shift of its centroid energy. 
Thus, some part of the fine structure of the experimental cross sections at lower energies is well described by the 
cPVC framework. The cPVC method based on a real NN effective interaction is found to successfully explain 
about 85% of the reaction cross section, through explicit channel-coupling effects. 
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Description of nucleon-nucleus (NA) elastic scattering 
based on the fundamental nucleon-nucleon (NN) interaction 
is one of the most challenging subjects of nuclear reaction 
studies, and is crucial for exploration of unstable nuclei, for 
which phenomenological optical potentials have not been es- 
tablished. The most essential quantity for this subject is the 
imaginary part W of the optical potential, which is responsi- 
ble for a loss of the incident flux due to the existence of non- 
elastic channels. Reliability of W can be judged by comparing 
the resulting reaction cross section or with experimental data. 

One of the most successful approaches to this goal is the 
folding model based on a complex effective NN interaction 
and a phenomenological or microscopic nuclear density. In 
Ref. Ql, for example, the folding model calculation with no 
free parameter is shown to well reproduce differential cross 
sections and spin observables for proton elastic scattering on 
12 C, le O, 40 Ca, 90 Zr, and 208 Pb at 65-200 MeV, as well as 
ctr of neutron on these targets at 20-800 MeV. The agree- 
ment is at almost the same level as that of the well-established 
Dirac phenomenology |00]. Many studies in this direction 
have successfully been done with the coordinate-space repre- 
sentation J3-01 an d the momentum-space representation 10] . 
Nowadays, the microscopic description of nucleus-nucleus 
(A A) scattering becomes a hot topic |9[. In all these mod- 
els, W comes from the imaginary part of the effective NN 
interaction, for which a Bruckner g matrix evaluated in infi- 
nite nuclear matter is often adopted. Despite the great success 
of this folding-model approach, it is quite obvious that the use 
of the g matrix is not feasible for low-energy scattering, be- 
cause individual energy-levels of the N + A system near the 
nucleon threshold strongly affect the scattering process. 

An alternative approach to the microscopic description of 
NA elastic scattering is many-body calculation of the N + A 
system including channel-couplings to the continuum states. 
In this approach, a NN effective interaction with no imagi- 
nary part is used, and W is generated through the coupling 
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to non-elastic channels that are explicitly taken into account. 
One of the most suitable models for this purpose will be the 
particle-vibration coupling (PVC) method, which describes 
collective vibrations and single-particle motion of individual 
nucleons simultaneously. In a recent paper 111011 . the micro- 
scopic continuum PVC (cPVC) method has been proposed 
and applied to studies on single-particle (sp) structures in 
40 Ca, 208 Pb, and 24 0. In general, the sp states of a nucleus 
A are observed in excitation energy spectra in the neighbor- 
ing Ail nuclei. It was shown in Ref. II 1 QTI that the cPVC 
method describes quite well the fragmentation of the sp hole- 
and particle-states as well as the shift of those centroid ener- 
gies, in good agreement with experimental data. These are 
distinguished features of sp states that are not taken into ac- 
count in the standard Hartree-Fock (HF) picture. 

The cPVC method is based on the self-consistent micro- 
scopic HF and the continuum random-phase-approximation 
(RPA) ifTTl [T2I1 with the Skyrme effective interaction. In this 
framework, the microscopic nucleon optical potential is char- 
acterized by the nucleon self-energy corresponding to specific 
energy E in the asymptotic region of the N + A system; E 
can be interpreted as the incident energy of the nucleon on the 
target nucleus A in the optical model picture. A great advan- 
tage of the cPVC method to other existing PVC models lfl3l - 
n~6fl is, as emphasized in Ref. IToh . the proper treatment of the 
continuum with imposing an explicit boundary condition on a 
nucleon in a continuum state. Therefore, it is quite promising 
that the cPVC method can be applied to studies on reaction ob- 
servables, for which boundary conditions play essential roles 
in general. 

In this Rapid Communication, we apply the cPVC method 
to the neutron scattering on le O below 30 MeV. We see how 
the cPVC method can describe the absolute values and en- 
ergy dependence of or in particular. Important aspects of the 
present study, in comparison with the preceding works iflrjl 
[1711 along the same line, are as follows. First, we treat the con- 
tinuum explicitly as mentioned above. Second, we include a 
very large model space (various phonon states up to 60 MeV) 
as described below. Third, we are interested in the structures, 
i.e., nontrivial energy dependence, of ctr due to the PVC, ex- 



2 



pected to appear at quite lower incident energies. 

In the cPVC framework, the scattering wave function of 
neutron v I'pvc( rcr : k) from A, with the relative coordinate r, 
the intrinsic coordinate a due to the spin degrees of freedom, 
and the relative wave number k in the asymptotic region, is 
described by the following Lippmann-Schwinger equation 
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(r<7, k) — 0f (ra, k) 

(ra, r'a'; E) 



x [v (rV) 5 (r' - r") 

+E(rV,//;£)]# (rV'.fc), 

(1) 

where 0f denotes the neutron free wave and v (r'cr') is 
the HF one-body mean-field potential. The PVC Green 
function and the corresponding self-energy are denoted by 
(ra, r'a'; E) and £ (r'a', r"a"; E), respectively. With 
the partial wave expansion, one may find that the transition 
matrix (T matrix) of the elastic scattering is given by 
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(2) 

where I (j) represents the orbital angular momentum (total 
single-particle spin) of neutron, and h\ (kr) and ji(kr) are, 
respectively, the spherical Hankel function with the outgoing 
asymptotics and the spherical Bessel function. Explicit ex- 
pressions of and £ are given by Eqs. (6) and (7) of 
Ref. ifloll . respectively. Note that we use E for the rela- 
tive energy between neutron and A, which is denoted by ui 
in Ref. iflQIl . Furthermore, we put (+) in the superscript of 
$pvc, G, and Gy, to explicitly represent that these func- 
tions satisfy the outgoing boundary condition. We solve a 
Dyson equation, Eq. (10) of Ref. [10], to obtain Gy , and 
hence T? vc . All other details can be found in Ref. ifToTl . 



The differential elastic cross section da/dtt is given by 
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Here, Pi and P( are, respectively, the Legendre polynomial 
and its derivative with respect to cos 9. The total cross section 
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The reaction cross section ctr is defined by 

a R (^)=a tot ( J B)- C 7 el ( J B). 



(8) 



As one may see from Eqs. ([TJ-ffJ]), the self-energy £ serves 
as a dynamical polarization potential. Therefore, the imag- 
inary part W of the optical potential comes from S, which 
takes into account the coupling within the Skyrme continuum- 
RPA response function and the HF Green function 1 10]. Note 
that we explicitly treat the nonlocality of E, whereas it was 
approximately localized in the previous work iflql . The HF 
T matrix can be obtained by 
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where means the HF Green function. 

In the present calculation, we adopt the Skyrme NN ef- 
fective interaction SkM* [18]. Note that, in contrast to in 
Ref. id 

, we do not introduce any other interactions to the cal- 
culation of scattering observables, i.e., a fully consistent treat- 
ment of the effective interaction is carried out. For the cPVC 
calculation, as in Ref. [10], the orbital angular momentum cut- 
off for the unoccupied continuum states is set at l cu t — 7h, and 
we include phonons associated with the multipolarities A w of 
2+, 3 _ , 4+, and 5", up to 60 MeV of the RPA excitation en- 
ergy. We choose the Fermi momentum fcp = 1.33 fm _1 for 
the residual force in the self-energy function of the PVC calcu- 
lation. The radial mesh size is Ar = 0.2 fm and the maximum 
value of r is set to 20 fm, to obtain the T matrix by Eq. 

First, we show in Fig. Q] the 77-dependence of a% t(E) of 
neutron scattering by le O; 77 is a parameter introduced in the 
evaluation of the Green function and the RPA response func- 
tion ilOll . and corresponds to the resolution scale of E. The 
dashed, dotted, solid, and dash-dotted lines show the results 
with 77 = 0.4, 0.2, 0.1, and 0.04 MeV, respectively. Consid- 
ering the experimental situation, E should be an order of eV, 
which cannot be achieved because of the computational lim- 
itation. As one sees from the figure, however, a tot (E) for 
E > 4.0 MeV converges at 77 = 0.1 MeV. Thus, in the fol- 
lowing discussion we use 77 — 0.1 MeV and restrict ourselves 
for E > 4.0 MeV. 

In Fig.Q] the solid line reproduces quite well the energy de- 
pendence of the experimental data [19] except for 10 < E < 
15 MeV. A remarkable feature is the reproduction of the fine 
structure of ctot (E) in part, which is the benefit of the PVC. 
It is well known that many peaks in the experimental atot(E) 
correspond to the compound resonance. It is quite obvious 



3 




10 15 
E [MeV] 



20 



25 



FIG 

16 



1: (Color online) Total cross section of neutron scattering by 
O, as a function of the incident energy. The dashed, dotted, solid, 
and dash-dotted lines show the results with r\ — 0.4, 0.2, 0.1, and 
0.04 MeV, respectively. Experimental data are taken from Ref. 11911 . 



that the present cPVC calculation cannot describe such peaks. 
Rather than that, the cPVC method is considered to well de- 
scribe the so-called doorway states [20]. Keeping this in mind, 
we see that the present calculation reasonably reproduces the 
energy dependence of crtot(^)- Note that we do not tune any 
adjustable parameters, although the result can, to some extent, 
depend on the Skyrme parameters adopted. 
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FIG. 2: (Color online) Total-elastic cross section of the neutron 
scattering by O. The solid (dashed) line shows the result of the 
cPVC method (HF calculation). Experimental data are taken from 
Ref. (H. 



Figure |2] shows the a c \(E) calculated by the cPVC method 
(solid line) and the HF calculation (dashed line). Results are 
plotted in the linear scale. The cPVC method reproduces well 
the experimental data 11911 . whereas the HF calculation over- 



shoots them by about 80% above 15 MeV. This shortcoming 
of the HF calculation is rather trivial because of the absence 
of the imaginary part of the optical potential; absorption ef- 
fects become more significant as E, the number of open chan- 
nels in fact, increases. It should be noted that the HF calcu- 
lation gives a c \{E) — atot(E) and hence o"r(E) = 0. An- 
other important finding is that the HF result has only one peak 
around 4.6 MeV, which turns out to be due to the £7/2 sp orbit. 
The cPVC method gives much more complicated shapes of 
a c \(E), as in Fig.Q] This result clearly shows the importance 
of the PVC in the neutron elastic scattering at low energies. 
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FIG. 3: (Color online) Reaction cross section of neutron for 16 0. 
The solid, dashed, dotted, and dash-dotted lines correspond to the 
calculation with A max = 5, 4, 3, and 2, respectively. Experimental 
data are taken from Ref. flf^l . 



Next we compare the result of <?r(E) (solid line) with the 
experimental data 11911 in Fig. [3] Again, we plot the result 
in the linear scale. Although the theoretical ctr(-E) slightly 
undershoots (overshoots) the data for 10 < E < 25 MeV 
(E < 6 MeV), some peak structure is reproduced well. On the 
average for 8 < E < 30 MeV, the discrepancy between the 
theoretical values and experimental data for ctr (-E) is about 
15%. It should be noted that the present cPVC calculation de- 
scribes the ctr(-E) only though channel-coupling effects, i.e., 
with no imaginary part of an effective interaction. It will be a 
remarkable achievement that about 85% of the <tr(-E) is suc- 
cessfully explained in this manner. It was shown in Ref. lfl7ll 
that more than half the <tr (E) for proton scattering on 58 Ni, 
48 Ca, and 90 Zr is due to a coupling with the deuteron (trans- 
fer) channel. In the present calculation, the deuteron chan- 
nel seems to play less important roles; it should be remarked 
that i) theoretical calculation shows good agreement with data 
for E > 25 MeV and ii) we still have undershooting around 
10 MeV, where the (n, d) reaction channel is closed. 

In Fig. [3] the dependence of a-R,(E) on the maximum mul- 
tipolarity A max is also shown. The dashed, dotted, and dash- 
dotted lines correspond to the calculation with A max = 4, 
3, and 2, respectively; the solid line represents the full cal- 
culation with A max — 5. One sees a good convergence at 
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Amax = 4. Another finding is the role of A = 3, which essen- 
tially generates the many peaks in <tr (E); A = 4 then gives a 
slight change in the shapes and positions of the peaks. 
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FIG. 5: (Color online) Angular distribution of the neutron elastic 
cross section by ie O at 28 MeV. The solid (dashed) line shows the 
result of the cPVC method (HF calculation). Experimental data are 
taken from Ref. d. 



FIG. 4: (Color online) Partial cross sections of cr to t. The solid 
(dashed) and dotted (dash-dotted) lines represent the results of the 
{7/2 and f 5 /2 orbits obtained by the cPVC method (HF calculation). 



To see how the PVC gives peaks in the cross section more 
clearly, we show in Fig.|4]the partial cross sections of otot(-E), 
otot;ij (E) in Eq. ©, for the f-orbits. The cPVC (HF) result 
for the m and f 5 / 2 orbits are shown by the solid (dashed) and 
dotted (dash-dotted) lines, respectively. The dashed line has 
a resonant peak around 4.6 MeV, which is fragmented by the 
PVC as shown by the solid line. Similarly, the PVC changes 
the dash-dotted line to the dotted line, generating nontrivial 
structures. Thus, the PVC completely changes the shapes of 
the partial cross sections. Fragmentation of the sp strength 
function is necessary to explain the energy dependence of ex- 
perimental data, as shown in Figs.Q~|-[3] 

Finally, we show the result of the differential cross section 
d<r/d£l in Fig. [5] The cPVC and HF results are denoted by the 
solid and dashed lines, respectively. The dashed line severely 
overshoots the data at backward angles. This issue remains for 
the solid line, although the agreement with the experimental 
data is significantly improved by including the PVC. This may 
indicate that the radial dependence of the optical potential is 
not correctly generated by the present framework; further in- 
vestigation on this is our future work. 

In this study, we applied the continuum particle-vibration 
coupling (cPVC) method to the calculation of scattering ob- 
servables, i.e., the total (<7tot)> total-elastic (<r e i), and reac- 
tion cross sections (or), of neutron on the le O target at 4- 
30 MeV. The cPVC method describes the single-particle mo- 
tion and collective vibrations simultaneously, with properly 
treating the boundary condition of a nucleon in continuum 
states. This feature makes the calculation of the transition 
matrix (T matrix) straightforward. We used the Skyrme in- 
teraction (SkM*) as an effective nucleon-nucleon (NN) in- 
teraction, and included phonon states of le O up to the 5~ 



state with excitation energy of 60 MeV. The T matrix is evalu- 
ated with the PVC Green function, which is obtained by solv- 
ing the Dyson equation. The couplings to various collective 
states are characterized by the self-energy S. Since the cPVC 
method employs a real NN interaction, the imaginary part 
of the optical potential purely comes from S. In the present 
framework, or is described as a loss of the incident flux to 
various channels that are explicitly taken into account. 

The results of the present cPVC calculation satisfactorily 
agree well with experimental data of otot an d <?eb For or, the 
cPVC method explains about 85% of the experimental data 
on average, which will be an important achievement of the 
mean-field type calculation for neutron scattering. Another 
remarkable feature of the cPVC result is the fragmentation 
of a single-particle resonant cross section. This results in a 
good correspondence with some of the peaks observed, prob- 
ably those due to the doorway states. On the other hand, the 
Hartree-Fock (HF) calculation was found to give a rather triv- 
ial shape of the cross section and severely overshoot a e \ for 
E > 10 MeV. The cPVC method describes the angular dis- 
tribution of the elastic cross section much better than the HF 
calculation. However, there still remains a discrepancy at scat- 
tering angles larger than 30°. This will be due to an incorrect 
radial dependence of the optical potential in the cPVC frame- 
work. 

In future, we will apply the cPVC method to other reaction 
systems, including proton scattering. Application to inelastic 
scattering and photo-induced reactions will also be important. 
Another important subject is the interaction dependence. The 
scattering observables can be used as new constraints on the 
parameters of the Skyrme interaction, determined to repro- 
duce nuclear bound-state properties. In addition to that, as 
discussed in Ref. II 1011 . further investigation on the treatment 
of Pauli's principle in the cPVC framework will be necessary. 
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